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Abstract: We introduce a universal method to optically induce mul- 
tiperiodic photonic complex superstructures bearing two-dimensional 
(2D) refractive index modulations over several centimeters of elongation. 
These superstructures result from the accomplished superposition of 2D 
fundamental periodic structures. To find the specific sets of fundamentals, 
we combine particular spatial frequencies of the respective Fourier series 
expansions, which enables us to use nondiffracting beams in the experiment 
showing periodic 2D intensity modulation in order to successively develop 
the desired multiperiodic structures. We demonstratively present the gener- 
ation of 2D photonic staircase, hexagonal wire mesh and ratchet structures, 
whose succeeded generation is confirmed by phase resolving methods 
using digital-holographic techniques to detect the induced refractive index 
pattern. 
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1. Introduction 

Inspired by numerous prototypes in nature, photonic structures have become a captivating re- 
search subject of great demand and luminous future perspective. Due to excellent theoretical 
as well as experimental results in the last decades, this field already educed many basic mod- 
ules for the future technology as to telecommunication and information industry, all-optical 
computing as well as displaying devices. 

The groundbreaking fundamental insight is the existence of structure- specific band struc- 
tures |[1 1, for instance providing complete band gaps that allow for particular Bloch modes to 
explicitly be reflected. In analogy to charge carrier propagation processes in condensed matter, 
many new effects of light propagation in refractive-index modulated photonic structures be- 
sides Bragg diffraction |2|, e.g. Bloch oscillation, Zener tunneling, and Anderson localization 
in the linear |j3j-|5J, as well as spatial discrete bright, dark, and vortex solitons in the nonlinear 
regime |[6]-[9|. 

Up to now, relatively simple systems were investigated, however the scientific interest as- 
pires at the investigation of more complex systems such as quasiperiodic |[TO|[TT|, curvilin- 
ear including Bessel and Mathieu lattices p2}fT4| where already solitary structures have been 
found, as well as randomized photonic structures p3| . One of the most interesting field of 
highly promising systems regarding spectacular propagation effects are multiperiodic photonic 
structures, as especially in these systems intriguing analogies of quantum mechanics such as 
Klein Tunneling or Zitterbewegung could be predicted theoretically as well as demonstrated 
experimentally |^16^J7 |, not disregarding distinct solitary solutions y_8J. 

In addition to direct-laserwritten |T9| or holographic-lithographically pQ| generated pho- 
tonic structures, photorefractive materials represent a highly appropriate system to carry out 
linear as well as nonlinear light propagation experiments [ f7]|21[22| . The main advantage here is 
the reconfigurable and dynamically manipulable induction of complex two-dimensional refrac- 
tive index modulations and thus, the optical induction of periodic, quasiperiodic and curvilinear 
two-dimensional [23 j as well as three-dimensional photonic structures could be presented |24|. 
Additionally, probing of the generated structures is feasible in one and the same setup and no 
further treatment of the sample is needed. 

The implementation of all the mentioned photonic structures were limited to one individ- 
ual structure-inducing wave field, that is generated by interferometry, diffraction on artificial 
holograms, or via wave field moulding per spatial light modulators. As it is advantageous to 
implement wide-elongated photonic structures offering a longer interplay length which offers 
an analogue of long distance propatation to temporal behavior, one is constricted to use writ- 
ing wave fields that are translation invariant in the direction of propagation, represented by 
the class of nondiffracting beams [25 J . Though, designing wave fields of multiple transverse 
periodic lengths by coherently summing up various nondiffracting wave fields each showing 
different structural sizes, the translation invariance of the intensity distribution of each wave is 
no longer conserved due to occuring longitudinal interference modulations. Anyhow, to imple- 
ment multiperiodic structures by optical induction using nondiffracting beams, one has to care 
for an incoherent superposition of several writing beams of various structural sizes, for instance 
by holographic multiplexing techniques that are well known from the field of holographic data 
storage ||26][27|. The technique of incremental multiplexing to modulate the refractive index of 
a photosensitive medium was presented earlier fSS", 29] where the former publication already 
presented the multiplexing method using 2D square lattices with different lattice periods as a 
proof of principle, and the latter demonstrates the generation of a one-dimensional (ID) pho- 
tonic ratchet structure. 



In this contribution, we propose an approach for the optical induction of two-dimensional 
multiperiodic photonic structures by incremental multiplexing techniques. The outline of the 
paper is as follows. After this first introductory section, Section |2] discusses theoretically the 
specification to create three representative multiperiodic structures. We respond to the proper 
incoherent combination of fundamental periodic lattices of different structural size to expand 
the desired superlattices. Recapitulatory, the mathematical considerations of each structure are 
summarized graphically. The experimental implementation of the developed series expansion is 
subject of Section |3j where the experimental setup is presented and details of the optical incre- 
mental multiplexing method are specified. In Section]?] the analysis of the generated photonic 
superstructures and the results are detailed. 

2. Multiperiodic structures and set of their fundamentals 

This Section addresses the derivations of the 2D superlattices in terms of Fourier series expan- 
sions, which can be later applied for the experimental part of optical induction. Therefore, we 
design the multiperiodic refractive index distribution by a stepwise and incremental induction of 
single index expansion terms using periodic nondiffracting beams (NDBs) with varying period- 
icity as the fundamentals. The usage of NDBs is essential, and thus we exclusively concentrate 
on this class of light waves as lattice-inducing wave fields and on their transverse intensity dis- 
tributions in particular. Although there is a wide spectrum of different intensity modulations 
among the four families of nondiffracting beams, we solely operate during the induction pro- 
cess with periodic discrete NDBs j30j[31| to preserve the periodic character of our photonic 
structures of inquiry. 

An intuitive approach to generate 2D structures with multiperiodic character is to use an 
optical Fourier series expansions via ID basis lattices in different orientations and varying 
periodic length. In this context, a 2D multiperiodic structure M is designed by combining ID 
lattice structures of equal modulation frequency resulting in a set of 2D fundaments of varying 
periodicity, which accounts for a series expansion of the desired structure. This procedure will 
for all of our examples end up in different sets of ten 2D basis lattice beams, where the series 
expansion is symbolically written as 

10 

M{r)=Y,amI{Kr). (1) 

m=l 

Here, am is the weighting coefficient of the mth order term and / describes a 2D basis intensity. 
Further, the vector r = (x, j)^ characterizes a point in 2D space and with j = 1 , 2 (or depend- 
ing on the rotational symmetry of the structure j = 1 , 2, 3, respectively) represents a modulation 
vector in the corresponding 2D Fourier space. In general, all following field distributions and 
their intensities are specified as 2D functions exclusively in a transverse plane. The thereto or- 
thogonal direction z can be identified with the direction of propagation. The field evolution in 
this direction as well as the temporal development are consequently neglected in the following 
considerations, as these behaviors can be examined by a multiplication of the field distribution 
of interest with the term exp(/(^||Z — cot)), where k\\ = {k^j^^ — k^)^^^ is the coordinate of the 3D 
wave vector in the direction of propagation and k describes the modulation in transverse 
direction. 

In the following, we present three representative multiperiodic structures with inherently 
different properties. Thereby, we concentrate on considerations of intensity distributions since 
with our technique exclusively the intensity is relevant for refractive index changes caused by 
the photorefractive effect. To achieve particular intensity distributions, we design corresponding 
light wave distributions showing nondiffracting properties as they can be very flexibly imple- 
mented experimentally |23|3H . Since the presented structures are only representative examples. 
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(a) (b) (c) 

Fig. 1. (Color online) Schematics describing a multiperiodic 2D staircase lattice; (a): sur- 
face plot of a simulated refractive index distribution generated by series expansion up to an 
order of ten, (b): presentation of the relevant spatial frequencies, (c) periodic square lattice 
as the basis structure of each expansion term. 



various combinations of rotational symmetry (even quasiperiodic) and expansion rules for arbi- 
trary periodic functions are implementable. Moreover, the combination of different expansion 
rules for unequal transverse directions are possible. 

2.1. Staircase superlattice 

The first structure we consider is a 2D staircase refractive index superlattice. This unique struc- 
ture offers steep slopes and flat plateaus, involving high spatial frequencies of nonnegligible 
significance. For that reason, a successful implementation of such a challenging structure proofs 
the high functionality and flexibility of the presented technique. 

Figure[TJa) represents graphically the desired 2D Fourier series expansion. In the illustrations 
of the relevant Fourier components sketched in Fig. [ijb), equally colored spots of a certain 
brightness represent the sine or cosine functions of one of the orientations. The first expansion 
orders are represented in gray boxes. Especially in these schematics, the multiperiodic character 
becomes obvious due to the higher order harmonics. Combinations of all contributions of one 
order m lead to the derived series expansions with 2D basis sets consisting of nondiffracting 
beam intensities. 

To derive the correct intensity distribution of each set of 2D fundamentals, we first analyze 
the expansion rule for a 2D rectangular function II2D with respect to ID fundamentals. To keep 
the intensity non-negative, we adjoin an offset term: 

10 Y 

^2D (r) = X [2 + sin (k» + sin (k^r) ] , (2) 

m=l 

where m determines the order of the expansion term. Thus, for every expansion term m two 
^m-modulations with orthogonal modulation vectors and k^ contribute to the total 2D dis- 
tribution Il2D- The length of the modulation vectors are given by the arbitrarily chosen lattice 
period of the first order g\ and the order m, hence: 

In 

k^ = (2m - 1) ^ e^- = km ej with j = 1,2, d = (1, 1), 62 = (1,-1). (3) 
^1 

In the next step we will combine the ID intensity contributions of a particular value of m. 
These terms reveal modulation vectors k^ which are transverse components of wave vectors 



with mutual equal projection length parallel to the direction of propatation and thus restrain 
any modulation of intensity in this direction which is characteristic for NDBs. The sum of the 
two sinusoidal intensity contributions for a fixed m in Eq. ^ can be rewritten as 



[2 + sin (201 ) + sin (202)] 



1 

'2 L 
1 



J(2(/)i-7r/2). 



-/(2(/)i-7r/2) , J(2(/)2-7r/2) _ 



-i(2(/)2-7r/2) 



(4) 
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where we introduced the phase function 0^- for a particular m with j = 1,2 using 20^^^ = k^r 
and dropping the index m. Now it becomes obvious that the wave field whose absolute square 
value appears in Eq. ^ equals the intensity of a field distribution consisting of four interfering 
plane waves. This field belongs to the four-fold discrete NDBs, whose transverse intensity 
distribution is presented in Fig.jTJc) and which is referred as ^4^1 (r) in Ref. |31|. In the case 
here, every plane wave carries an additional phase shift of ±;r/4, where plane waves of opposite 
spatial frequency have opposite sign, resulting in a real space translation of the structure in 
transverse direction (cf. Ref. |31 1). 
Hence, we can transcribe Eq. ^ to 



10 ^ 

2(2;;^ 



1^4,1 (r,/^m)r 



(5) 



and receive a series expansion of a multiperiodic 2D staircase structure incorporating transverse 
intensities of a nondiffracting wave field as the fundamental structures. 



2.2. Hexagonal wire mesh superlattice 

In this Subsection we derive an expression for the set of basis lattices creating a hexagonal 
multiperiodic structure. Due to their hexagonal property such systems are of particular impor- 
tance and coveted subject of many investigations for instance in graphene physics |32|. Basi- 
cally, a more general interest in hexagonal structures arises from its non- separable character 
which implies the impossibility to generate a hexagonal lattice by a sum of two orthogonal 



one-dimensional lattices. In analogy to the considerations of Subsection 2.1 we start with the 
sum of ID cosine lattices, resembling the hexagonal superstructure illustrated in Fig.[2|a): 



10 ^ 

^2D (r) = X - [3 + cos (k» + cos (k^r) + cos (k^r) ] . (6) 

m=l ^ 

In this case, we need a hexagonal basis set of modulation vectors, which becomes apparent in 
the illustration of the relevant Fourier components in Fig.|2jb): 

In 

kj^ = m — ej = km with j =1,2,3 

^1 

and 

(7) 

ei = (cos(;r/6),sin(;r/6)), ^ ^ 

e2 = (0,-l), 

63 = (— cos(;r/6),sin(;r/6)). 
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(a) (b) (c) 

Fig. 2. (Color online) Schematics describing a multiperiodic 2D wire mesh lattice; (a): 
surface plot of a simulated refractive index distribution generated by series expansion up to 
an order of ten, (b): presentation of the relevant spatial frequencies, (c) periodic hexagonal 
lattice as the basis structure of each expansion term. 



Again we use the abbreviation (j)j with in this case 0^^^ = k^r and j = 1,2,3 valid for all 
orders m. The term in square brackets appearing in Eq. ([6]) and representing the unweighted 
intensity of every expansion order m can be rewritten to 

[3 + cos (01 ) + cos (02) + COS (03)] = ^ [6 + ^^'^^ + ^"^'^^ + + + + ^"^'^'] . (8) 

Now we have to find a 2D basis for the expansion, for instance given by the intensity of a 
hexagonal nondiffracting beam, which consists of three plane waves conforming the spatial 
frequency conditions of a nondiffracting beam (cf. Ref. |[3Tj): 



|^3,o(r,/:m)|' 



(9) 



which can be identified with Eq. ([5]) using 0i = (pi — (P2, 02 = 9i — 93, 03 = <P2 — 93- Now, the 
arguments of the three plane waves in Eq. ^ are summations of the formerly introduced phase 
functions 0^. 

Thus, the rule for the development of a multiperiodic 2D wire mesh function in hexagonal 
symmetry reads as 



10 



(10) 



and is equivalent to Eq. ([6]). 

Besides the weighting of the fundamental structures, also the choice of the rotational sym- 
metry brings out a parameter to vary the resulting multiperiodic structure. Thus the six-fold 
periodic structure presented in Fig. |2jc) is the basement for the hexagonal multiperiodic wire 
mesh structure presented in Fig.[2|a) which conserves the rotational symmetry. 

2.3. Ratchet sup erlattice 

In general, functional structures such as 2D ratchets have various intriguing applications in 
research fields such as biological microscopic machinery and Bose-Einstein condensation [ 33 1, 
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(a) (b) (c) 

Fig. 3. (Color online) Schematics describing a multiperiodic 2D ratchet lattice; (a): surface 
plot of a simulated refractive index distribution generated by series expansion up to an 
order of ten, (b): presentation of the relevant spatial frequencies, (c) periodic square lattice 
as the basis structure of each expansion term. 



as well as quantum mechanics in solid-state or atomic physics systems 1 34] due to their ability 
for long-distance mass transport caused by a low-amplitude potential. 

In analogy to the considerations regarding the two foregoing structures, we are going to 
develop the rule for the optical induction of a 2D ratchet structure. Therefore, we again start 
with a Fourier series expansion under the usage of a basis set of ID fundamental lattices: 

10 ^ 

A2Z)(r) = ^ - [2 + sin (k» +sin (4r)] . (11) 

m=l ^ 

Figure|3ja) depicts the multiperiodic lattice in terms of a simulated refractive index distribution. 
Analog to Fig. [TJb) the corresponding series expansion respects specific spatial frequencies 
which are depicted in[3jb) for the first three orders of expansion. 

The set of modulation vectors assembles for each order m from two mutual orthogonal vec- 
tors and are chosen to be 

27r 

ki = m^ej = ej with j = 1,2, d = (1,1), 62 = (1, -1). (12) 
^1 

With the previously introduced phase function 0^ = kjr/2, where 7 = 1,2, and under consider- 
ation of Eq. (|4]), every expansion term of A2D can be expressed through the earlier introduced 
wave field ^4j(r,^^), which again is depicted in Fig. [3];c). Hence, the series expansion for 
a multiperiodic 2D ratchet intensity distribution A2D composed of a 2D set of fundamental 
nondiffracting beam intensities is 

10 I 

A2D(r)= 1^ ;^|^4,i(r,U^ (13) 

where exclusively the weighting factor 1 / (2m) as well as the relation between the modulation 
vectors of every order m determines the difference between Eq. ^ and Eq. ( 13 ). 



3. Experimental implementation 



In the previous Section, we found the rules for the generation of multiperiodic structures deriv- 
ing basic structures which resemble to intensities of nondiffracting beams. With these findings. 
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Fig. 4. (Color online) Schematic of experimental setup to optically induce and analyze 
spatial refractive index modulations of photorefractive crystals. Green colored light path 
represents the lattice inducing part, the probing part is colored in red. A/PSLM: ampli- 
tude/phase spatial light modulator, BS: beam splitter, D/FSC: direct/Fourier space camera, 
L: lens, A/2: half- wave retardation plate, M: mirror, MO: microscope objective, P: polar- 
izer, PH: pinhole, SBN: cerium-doped strontium barium niobate. Camera picture presents 
a typical recorded intensity pattern, magnification shows details of interference fringes. 



we are able to experimentally implement corresponding photonic structures, where the induc- 
tion procedure orientates on the theoretical series expansion mentioned in equations ([5]), ([TO]) 
and ([13]). Here, the usage of 2D fundamentals bears an essential advantage regarding the tem- 
poral effort in comparison to ID basic structures since for a series expansion of up to ten terms 
30 ID lattice structures are needed, which differ among each other in modulation frequency 
(i.e. lattice period) and orientation, as well. Using 2D fundamentals, we only need 10 different 
multiplexing steps. 

The setup scheme for the experimental implementation of the optical induction of multiperi- 
odic photonic structures is illustrated in Fig.|4] An expanded continuous-wave solid-state laser 
beam with wavelength A = 532 nm is phase modulated by a computer controlled spatial light 
modulator (PSLM). The phase pattern given to the PSLM and the displayed pattern of the fol- 
lowing amplitude modulator for spatial frequency filtering (ASLM) are chosen to implement 
particular nondiffracting beams according to certain expansion terms. The phase distribution 
includes both amplitude as well as phase information about the desired wave field which is 
numerically calculated beforehand. Besides phase modulation, amplitude modulation can be 
achieved by an intensity weighted blazed grating to diffract the relevant light into a certain 
diffraction order |35|. Experimentally, the illumination time tm is easily adaptable via a short 
switching time of the PSLM. 

By means of further optical elements, a desired wave field - in our particular case a non- 
diffracting beam - is imaged into a volume of interest, where a photorefractive crystal is placed 
in order to modulate its refractive index. We employed a nonlinear cerium-doped strontium 
barium niobate (SBN) crystal that, externally biased, translates the intensity modulations into 
refractive index changes via the photorefractive effect, as earlier described in |23|. For the 
whole set of induction patterns, the external voltage is about 800 V over a distance of 5 mm 
where the induction beam propagates along a 5 mm orientation of the crystal. 

As a generalized approach, the temporal development of the refractive index change in a 
biased SBN crystal caused by an illumination pattern behaves exponentially [27 J , fades to sat- 



uration with the time constant Ty^, and also gets exponentially erased (t^) due to illumination 
with non-identical intensity distributions. This can be summarized by 



with A^o = ^n{t = 0) and ^Usat as the saturation value of the refractive index change. 

In all multiperiodic expansions introduced in the previous section, the intensity belonging to 
a particular expansion order m is weightened by a factor depending on m. This weighting factor 
is realized in the experiment by adapting the illumination time tm for each lattice-inducing 
beam. Thereby, rather than controlling the external voltage or intensity of each fundamental 
structure, we regulate tm to weight a specific induction term. Hence, we determine the net 
illumination duration for one sequence of illumination to 4eq = Ll=i^m = 10 s and execute 
about 20 multiplexing sequences. 

The optical path of the index-modulating part in Fig.|4]is colored in green, whereas the red- 
colored structure-analysing part is subject of the next section. 

4. Analysis of photonic structures 

In order to analyze the optically induced photonic structures, we apply digital-holographic tech- 
niques in terms of superimposing a probe beam propagating through the generated structure 
with a reference beam of high spatial frequency. The detected phase- sensitive interference pat- 
tern allows to determine the whole wave field information using a common CMOS camera 
rather than to merely detect the intensity distribution. In this way the phase retardation of a 
probe beam propagating through a refractive index modulated crystal gives 2D resolved infor- 
mation about the refractive index structure and depth |36|. 

To extract the phase information of the probe beam from the recorded interference pattern, 
one has to treat the data in analogy to an amplitude demultiplexing operation in electric signal 
processing. According to this, bandpass filtering in the vicinity of the carrier frequency followed 
by a frequency shift to zero provides the complex field information of the signal which is equal 
to the probe wave field and includes the desired phase information p7| . In order to bandpass 
filter the signal, we multiply the spatial frequency spectrum with a radial symmetric Hann 
window, whose center is the carrier frequency and its radius equals the half distance between 
carrier and zero frequency. In addition to a real space camera, we implemented a further camera 
to detect the intensity of the probe beam in the far field by usage of the Fourier-transforming 
attribute of an additional lens. 

Figure|5]presents the results of the analysis of the induced photonic staircase [(a)-(c)], hexag- 
onal wire mesh [(d)-(f)] as well as ratchet structures [(g)-(i)]. Intensity as well as phase dis- 
tributions of the probe plane wave in real space combined with the intensity distribution in the 
far field are illustrated in the lefthand and middle column, respectively. In order to determine 
the refractive index change of a particular area, we first identify the phase distribution of the 
homogeneous medium before the induction process was started. This distribution serves as the 
reference phase and is subtracted from the measured phase distribution of the induced struc- 
ture. In this manner, exclusively the refractive index change A^ caused by the optical induction 
process can be extracted by carrying out An = A(pA / (2;rJ), where d is the optical path length 
of the probe beam through the medium. 



The righthand column of Fig. [Slvisualizes in a surface plot the refractive index landscape of 

an area of approximately 0.16 mnr. In addition to these illustrations, the induced structures are 
presented as animated surface plots supporting an enhanced spatial impression of all refractive 
index modulations (cf. Media 1, Media 2, Media 3). 
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Fig. 5. (Color, Multimedia online) Analysis results of probing experiments using propa- 
gating plane wave of perpendicular incidence to crystal surface to characterize generated 
2D photonic superlattices: (a)-(c) staircase structure, (d)-(f) wire mesh structure, (g)-(i) 
ratchet structure. Direction of c axis is illustrated in (g) and (i), respectively. Scale bar for 
all real space pictures is appended in (h). Left column present intensity, middle column 
phase distribution of the plane probe wave in real space. Insets of left column pictures il- 
lustrate far field diffraction spectrum. Refractive index landscapes are depicted as surface 
plots in right column. Animated plots for photonic staircase (Media 1), wire mesh (Media 
2) as well as ratchet structure (Media 3) are available online. 



Both the intensity and phase distribution of all three propagation measurements reveal the 
character of the respective structures introduced in Section |2] The staircase superlattice pre- 
sented in Fig. [5ja)-Fig. |5jc) is in great agreement with the simulations (cf. Fig. [T]) showing 
three kinds of plateaus of distinct level where the intermediate level appears most [which be- 
comes apparent in the surface plot of Fig.[5jc)]. The structure is thoroughly four-fold as the 
real-space results and the far field intensity distribution indicate, as well. According to this, the 
diffraction orders given in the inset of Fig. |5ja) resemble a discrete diffraction pattern multi- 
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Fig. 6. (Color, Multimedia online) Time development of the rectangularly modulated re- 
fractive index structure. The detailed illumination times are (a) ti = 20 s, (b) t2 = 80 s, 
(c) ^3 = 140 s, (d) t4 = 200 s. Respective insets present corresponding intensity distribu- 
tions in far field. The temporal development of the staircase structure can be found online 
as a video file (Media 4). 



plied by a 2D four-fold sine function. We note that this specific rotational symmetry particularly 
results from the diamond orientation of the fundamental structures and is not the ordinary case 
for photorefractive media with drift dominated charge carrier transport processes causing an 
orientation anisotropy p8| . 

In contrast to the staircase superlattice, such a anisotropy becomes apparent for the hexagonal 
wire mesh structure where intensity and phase distribution [cf. Fig. [5jd)-Fig. [5jf)] display a 
stronger refractive index modulation along the crystal's symmetry axis (c axis) compared to 
any other direction. Thus, vertical lines perpendicular to the c axis are more pronounced than 
the ones that are inclined by 30° to the c axis. The far field picture presented in the inset of 
Fig. ^d) and the surface plot in Fig. [Sj^f) emphasize the influence of the crystal's orientation 
anisotropy, additionally. 

With the ratchet superlattice presented in Fig. [5jg)-Fig. |5ji), we introduce a further 
anisotropy-unaffected structure besides the staircase superlattice, though this structure is not 
rotational symmetric as the straightering character is inherent for the ratchet structure. In par- 
ticular, the far field intensity illustrated in the inset of Fig. |5jg) manifests a preferred direction 
of light diffraction (spectral quarter above center). Moreover, in contrast to the ID ratchet struc- 
ture presented in |29|, the 2D ratchet superlattice can be oriented in four different ways without 
disturbance by the anisotropy and can further be used to direct light parallelly as well as antipar- 
allelly or perpendicularly up as well as down to the c axis, respectively. Altogether, especially 
the ratchet structure is preferred to be applied as a light director and marks a throughout relevant 
system for propagation experiments in the linear and in the nonlinear regime, as well. 

In the setup presented in Fig.|4]we use two distinctive wavelengths - 532 nm for the induction 
of the refractive index modulation (frequency doubled solid state Nd: YAG cw laser) and 633 nm 
as the probe beam wavelength (He-Ne cw laser). The advantage of using red laser light in our 
probing experiments is the smaller sensitivity of SBN in this wavelength regime compared to 
wavelengths in the green regime around 532 nm |28J . 

By using a separate probe laser, we are able to timely-resolved detect the refractive index 
modulation, which facilitates the development of the index change due to optical induction. 
Figures [6ja)-[6jd) present the developing index structure of a fixed area for four different points 
in time ti = 20 s, t2 = 80 s, = 140 s, and t4 = 200 s. The formation of the staircase pho- 
tonic structure becomes obvious already after = 140 s, and also the far field intensity reveals 
characteristic details of the desired structure after several seconds of illumination. 



5. Conclusion 

In conclusion, we introduced a ubiquitous method to generate 2D photonic superstructures by 
use of 2D nondiffracting beams. These writing beams serve as fundamentals of an optical se- 
ries expansion to generate multiperiodic photonic structures. In such a manner, e.g. anisotropic, 
graphene-like, as well as discontinuous functional structures are realizable. Although we proved 
a successful induction of three representative complex superlattices in a photorefractive SBN 
crystal, the technique is easily transferable to all photosensitive media, enabling the develop- 
ment of photonic structures of long expansion in transverse as well as in longitudinal direction. 
Also the implementation of multiperiodic optical atom or particle traps is one prospective appli- 
cation of this technique. In general, our method offers a highly-functional tool for the generation 
of two-dimensional multiperiodic systems, providing the basement of the investigation of so far 
unexplored fascinating linear as well as nonlinear light propagation effects. 
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